setwd("D:/Rev2Plot/")
rt<-read.table("AllPara.dat",header=TRUE)
mu<-2.33
m_H<-1.67e-24
kB<-1.38e-16
m_18<-(12+18)*mu
kms<-1e-5
pc=3.086e18
msun<-1.99e33

options(digits=2)
i=1
while(i<=length(rt$Name))
{
    if(rt$FW18[i]!=0)
    {
      FWHM18<-rt$FW18[i]
	    Tex<-rt$Tex[i]	
	    R<-rt$R[i]
	    n<-rt$n[i]
	    sTH<-rt$sTH[i]
      rt$sTH[i]<-sqrt(kB*Tex/m_H/mu)*kms   
      
    	sNT<-FWHM18/sqrt(8/log(2))
    	sNT<-sNT^2-kB*Tex/m_18*(kms^2)
    	
    	if(sNT>0)
    	{
    		sNT<-sqrt(sNT)
    	}
    	else
    	{
    		sNT<-0
    	}
    rt$sNT[i]<-sNT
    	
    	sTD<-sqrt(3*(sTH^2+sNT^2))
    rt$sTD[i]<-sTD
    
    Mvir<-2.1e2*R*(FWHM18^2)
    rt$Mvir[i]<-signif(Mvir,2)
    
    Cs_eff<-sqrt(sNT^2+rt$sTH[i]^2)/kms
    Teff<-Cs_eff^2*mu*m_H/kB
    Mj<-(Teff/10)^1.5*(n*1e3/1e4)^(-0.5)
    rt$Mj[i]<-signif(Mj,2)
    }
    
    else
    {
      FWHM13<-rt$FW13[i]
	    Tex<-rt$Tex[i]	
	    R<-rt$R[i]
	    n<-rt$n[i]
			sNT<-rt$sNT[i]
    sTH<-rt$sTH[i]
    Mvir<-2.1e2*R*(FWHM13^2)
    rt$Mvir[i]<-signif(Mvir,2)
    
    Cs_eff<-sqrt(sNT^2+rt$sTH[i]^2)/kms
    Teff<-Cs_eff^2*mu*m_H/kB
    Mj<-(Teff/10)^1.5*(n*1e3/1e4)^(-0.5)
    rt$Mj[i]<-signif(Mj,2)
    }
    
    NH2<-rt$NH2[i]*1e21
      R<-rt$R[i]*pc
    Mlte<-NH2*R^2*pi*m_H*mu/msun
    rt$Mlte[i]<-signif(Mlte,2)
    
    i<-i+1
}   

setwd("D:/Rev2Plot/")                    
rt<-read.table("AllPara_New.dat",header=TRUE) 




region="CMC"

i=1
while(i<=length(rt$Name))
{	
	j=2
	if(rt$Region[i]!=region)
	{
		while(j<ncol(rt))
		{
			rt[i,j]<-NA
			j<-j+1
		}
	}
	i<-i+1
}

